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SUMMARY AND CONCLUSION * _ LV- - 

This paper investigates the mathematical relationships between 

the shape parameter and estimates of reliability and a life limit 
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bound and the allowable life limit lower bound (often called a 
tolerance limit) have unique global minimums over a range of Jr. 
Hence lower bound solutions can be obtained without assuming or 
estimating j3T "The existence and uniqueness of these lower bounds 
are proven _ in the Appendix. Some real data examples are given to 
show how these lower bounds can be easily established and to 
demonstrate their practicality. The method developed in the paper 
has proven to be extremely useful when using the Weibull 
distribution in analysis of no-failure or few-failures data. The 
results are applicable not only in the aerospace industry but 
anywhere that system reliabilities are high. 

INTRODUCTION 

The two parameter Weibull distribution is widely used for 
reliability and life limit estimations. However, there is often 
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insufficient information, e.g., failure data, to accurately 
estimate the shape parameter p of the distribution. Traditionally, 
either a value of -P is assumed and then used to establish the scale 
parameter a and estimate a reliability lower bound, or a 
nonparametric method is invoked to estimate a reliability lower 
bound. In the former situation the assumption of p is often 
controversial and difficult to justify. In the latter case, the 
reliability bound estimated is usually too conservative. This paper 
presents an alternate method which attempts to overcome the 
drawbacks in the traditional methods. 

NOTATION AND ASSUMPTION 

Notation: 

Cdf Cumulative distribution function 

p Weibull distribution shape parameter (slope value) 

a Weibull distribution scale parameter (characteristic 

value) 

p Estimate of any parameter, p 

n Total number of test units 

r Number of failed units 

ti Operating time (seconds, cycles) 

t m max(ti, t 2 , • • • , t n ) 

Zi tf/tm, normalized time 

T Specific operating time or life limit lower bound 

R Reliability for a given life limit 

T c Conditional life limit lower bound 

R c Conditional reliability 
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MTBF 


Mean time between failure 


r Confidence level 

Ln Natural logarithmn 

Assumption: 

For the hardware units, the distribution of operating time to 
failure is the 2-parameter Weibull. 

MATHEMATICAL RESULTS 

The 2-parameter Weibull Cdf is: 

F(t) = 1 - EXP ( , t>0 (1) 

Where t is a variable generally given in terms of time or cycles, p 
is the shape or slope parameter of the distribution and a is the 
scale or characteristic life parameter. Assume n units are tested. 
Let ti,t 2 ,...,t r be failure data and t r +i,t r + 2 , . • . ,t n 
be censored data (either failure censoring or time censoring) . 

A lower bound on a results from the fact that if t is 
distributed as Weibull {(3 , a) then tP has an exponential 
distribution with MTBF parameter aP . Thus given a value of p, a 
100t percent confidence lower bound on a (Refs. 1-4) is given by: 


; = ( f . ,1 /p 

C(T) 


where c(r) is a positive constant which depends on the confidence 
level t, the number of failures r and the censoring scheme. For 


example, for type II censoring, c(r) = •5*X J 7 - / 2r' f° r tYP e 1 
censoring with replacement, c (r) = . 5*X : T/ 2r+2 (Ref. 4), where 
X 2 T t k is the lOOr percentile of the Chi-square distribution with 
degrees of freedom k. 
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With high reliability systems a value of p is frequently hard to 
come by. Either there are no failures in which case p cannot be 
estimated or there are few failures in which case the estimate of p 
is subject to considerable uncertainty, and no assumption about p 
based on similar failure modes or similar hardware is available. 
Fortunately, it turns out that under certain conditions, which are 
not so restrictive, reliability for a given life or life limit for a 
given reliability experiences a minimum value as a function of p. In 
other words there is a "worst case (3", say Pq, such that 
reliability or life limit attains a strict minimum when P=Pq. The 
equations and conditions are shown below. 

Reliability lower bound: 

When the unit is tested to time T, it is easily seen that the 
100 r percent confidence lower bound on reliability is 

- - rpR T B 

R s R(jS) = EXP ( - -pr) = EXP( - c ( r ) stjL ff ) (3) 

Theorem 1 : Under the condition: (tit 2 • • • tn) -*-/ n < T < 

max(ti, t 2 , . . • , t n ) , there exists a unique p q>0 such that 

R (P) > R(Po) , for all p > 0 and p Pq; 
and R (P) is a monotone increasing function of p when P > P q, 
and R(/3) is a monotone decreasing function of p when p < Pq. 

For all theorems presented here, when the conditions are not met, 
the function in question is either monotonically increasing or 
monotonically decreasing over the range of p, in either case a worst 
case p does not exist. 

Allowable life limit lower bound: 

Assuming the reliability is a given constant R, we solve for T 
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from (3) , then a 100 r percent confidence lower bound on allowable 
life limit (tolerance limit) is 

T = T(/3) = ( - Ln(R)- |^ jr y - ) V/? (4) 

Theorem 2 : Under the condition: k < c(r)/ (-Ln(R) ) < n, where k = 
number of ti's with values = max (ti , t 2 , . . . , t n ) , 
there exists a unique /3o>0 such that: 

T(/3) > T (/?o) / for all jS > 0 and 0 + Po't 
and T (/?) is a monotone increasing function of f3 when /3 > /3o, 
and T(/3) is a monotone decreasing function of /? when /? < /Sq. 


Rocket engine reliability is more commonly thought of in terms 
of mission reliability rather than life time reliability. Thus 
there is interest in the reliability of the next to second 
mission given that the unit has already accumulated T-to seconds. 
The equations and conditions for the conditional reliability and 
life limit are given below. 

Conditional reliability lower bound: 

We define the conditional reliability R c as the probability 
that the unit will not fail during the next to seconds, given 
that the unit has survived the first T-to seconds. Then R c is 
computed by 


RC = RXP( - 


s? >■ <5 > 

Substituting equation (2) into equation (5), we get the lOOr percent 
confidence lower bound on R c : 


Rc = Rc(0) = EXP( - c(t) — } (6) 

Theorem 3 : Under the conditions: 0<to<T<max(ti, t 2 , . . • , t n ) 
and Sti (l+Ln(T/ti) ) <0 there exists a @o>l such that 
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' to/3nT^® - ^ 

R C (^) > EXP(-c (O — Q £ g" ), for all jS>1. (7) 

Et i 

Conditional operating life limit lower bound; 

We frequently need to assess a lower bound for a conditional 
operating life limit T c which is defined by (5) for a given 
conditional reliability R c . While T c can not be explicitly 
solved for from (5) , it can be treated as an implicit function of 
other variables. 

Theorem 4 . Under the conditions: 

R c < EXP ( - c(r)|^r-) and Eti*Ln(^j) > 0, (8) 

there exists a Po > 1 such that 

T C (0) ^ ( -Ln (R c ) c ^t ) ^Q t Q K (/?0 " 1) , for all p> 1. (9) 

Note that Theorems 3 and 4 have been restricted to /3>1 (wear 
out type failures) . Similar work for /3<1 (infant mortality type 
failures and constant failure rate type failures) will be the 
subject of a future paper. 

The proofs of the Theorems 1 to 4 are given in the Appendix 
and a computer program has been written to solve for the /3o's 
and the minimums . 

ILLUSTRATIVE EXAMPLES 

Example 1. Of the 59 units of a particular component that were 
tested only one failed. The operating times (seconds) of the 59 
units are: 14176 (failure), 28587, 20359, 18248, 17256, 16769, 
13890, 13357 13118, 13118, 13061, 12893, 12640, 11878, 11858, 

10378, 8172, 7918, 7379, 7326, 6700, 6442, 6263, 6168, 5809, 5383, 
5305,4788, 4788, 4784, 4748, 4701, 4453, 4258, 3696, 3145, 3138, 
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3059, 2867, 2837, 2777, 2438, 2300, 2178, 2027, 1807, 1739, 1567, 

1384, 1357, 1303, 1161, 1062, 1062, 1042, 903, 819, 299, 299. 

We wish to use the Weibull distribution to assess the 
reliability lower bound at the 50% confidence level for a new unit 
operated for 10,000 seconds. 

Since we have just one failure, we are not able to obtain a 
good estimate of /3 . Using Theorem 1 with: 

r = .50; c(r) = X 2 T/ 4=3.357 and T = 10,000 seconds, 

we have R(/3) = EXP( - . 5*3 . 357 ^|^^ ) 

From Theorem 1, we know there exists a unique Po>0 (Figure 1) such 
such that 

R(j8) > R(/3o) for all Po>0. 

The computer program computes: Pq = 1.083 and R(/3o) = -958. 

Therefore we may conclude with 50% confidence that the reliability 
of a unit operated for 10,000 seconds is at least .958. 

Example 2. Sixteen units of a particular component experienced the 
following cycles with no failure: 35, 31, 26, 21, 40, 30, 38, 27, 

23, 38, 18, 19, 12, 14, 13, 9. Given these data, what is the 
minimum number of tests which can be run without failure at .95 
reliability and 50% confidence level? 

We have: R=.95, r=.50 and c ( r ) =. 5*X 2 T f 2=1 • 386 . Using Theorem 
2 (see Figure 2) , we get: 

/3 o = 1-529 and T(/?o) = 28.6 starts 

So we conclude that at the 50% confidence level and 0.95 reliability, 
the component could be tested for at least 28 cycles without failure. 
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= 59 
' = 1 

' = 10 , 000 sec. 
=- .50 


1.083 / 

.e of Theorem 1 


Example 3. Using the data in Example 1 compute the 90% confidence 
lower bound on reliability for 520 seconds of operation given that 
the unit had already operated for 5000 seconds without failing. 

We have T=5, 000+520=5520 seconds, to=520 seconds, r=.90 and 
c ( r ) = . 5*X 2 T t 4 = 7 . 779 . Using Theorem 3 (see Figure 3), we get: 

Po = 1.321 and equation (7) gives 

R c (/3) > EXP ( -.5*7.779 520 * 1,321 *^ 5 ^^ 321 ~ 1 ) = .9945, for all 0>1, 

Eti 

Thus we conclude, with 90% confidence, that the conditional 
reliability of a unit operating for 520 seconds given that it had 
already operated for 5000 seconds is at least .9945. 

Example 4. Again using Example 1 data compute the life limit such 
that a 520 second operation reliability of at least .99 is 
guaranteed at the 95% confidence level. 

We have to=520 seconds, r = . 95, c (r ) = . 5*X J j 4 = 9 .4 88 and 
Rc=.99. Using Theorem 4 (see Figure 4), we get: 


Po = 2.188 and equation (9) gives 


T c (/3)>(-Ln(.99) 


gtj 2. 188 , 

.5*9. 488*2. 188*520 ; 


1/(2.188-1) 
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=9,486 seconds, for all (3>1. 


So we see, in order to meet 0.99 reliability lower bound at 95% 
confidence level, the conditional life limit is at least 9,486 
seconds . 




APPENDIX 


1. Proof of Theorem 1: 

T/3 

From (3), we get: R(/3) = EXP( -c ( r ) 

It is obvious that finding a minimum for R(j8) is equivalent to 
finding a minimum for g(/3) = E(ti/T)0. 

The first derivative and the second derivative of g(/?) are 

g' U3) = E Ln(| i )*(| i ) /? (1.1) 

and g"(/3) = E (Ln^) 2 * (i^) ^ (1-2) 

From (1.2) , and the given condition: 

(tit 2 . . . t n ) 1/n < T < max(ti,t 2 , - • • ,t n ) , (1.3) 

we see g"(/?)>9. So g'(/3) is a strictly monotone increasing function 
of /3. Also because of (1.3), we have 
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Jim o g'(0) = ELn(^) = Ln(^ . .^) < 0 

and ^in| +oo g 1 ()3) =+oo>o . Therefore there exists a unique (3o>0 such 
that : 

g'^o)- 0 and g (P) attains a minimum at f3 o; 

and g(/3) is a monotone increasing function of /3 when (3>Po, 
and g(/3) is a monotone decreasing function of /3 when P</3q. 

Transferring the result from g(/3) back to R(/?) / we conclude 
the Theorem 1 . 

2. Proof of Theorem 2: 

From (4) we get T(/3) = ( -Ln(R) 

Using the notations tm=max (ti , t 2 , • • . ,t n ) , zi=ti/t m , 
i=l,2,...n; and C=-Ln (R) /c ( r ) . We have 


%£>- (crzi/3) 1 ^ = exp ( totOttoas i a , 
tm h 


( 2 . 1 ) 


Define h(/3) s ( 


, Ln fC) +Ln (Ez i ) 


). We have = EXP(h(/3)). 

r m 


The first derivative of h(/3) is 


( 2 . 2 ) 


h'(j3) = ( -LnC + ~ Ln(E 2i /3) ) / 0 2 

Defining hi(/3) = -LnC + /3 SLn ^ “ Ln(Ezi/3) , 

we have h' (/?) = h]_(/3)//3 2 . The first derivative of hi(/?) is 

k i / o \ _ SLn ( z i Uil , r (2 (Lnz i ) 2 z-i ft) (Ez ;j ft) - (Ln ( z j ) z jft ) 2 
hl (/3) - Ezift + ^ (Ezift) 2 

ELn ( z j ) z 2 ft 
SziP 

^ ( E ( Ln z i ) 2 zift) jZzjll -(SLn(z-j)zift) 2 
- ft 

Using the Cauchy Inequality (Ref. 6, Page 11), we have hi'(/3)>0. 
But by the given condition k<c(r)/ (-Ln(R) ) <n, where k=number of 


(2.3 
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ti's with values=t m , the equality does not hold. So 
hi' (P)>0. Hence hi(/3) is a strictly monotone increasing 
function of p. 

Using the condition k<c (r ) / (-Ln (R) ) <n again, we have 
^li^ hi(/3) = -LnC + 0 - Ln (n) = -Ln(Cn) < 0 
and 0 lim hi(/3) = —LnC + 0 - Ln(k) = — Ln(Ck) > 0. 

p — ►-poo 

Therefore, there exists a unique Po>0 such that hi(/3o)=0/ 
and hi(/3)<0 when P<Po, 

and hi(/3)>0 when P>Pq. 

Noticing hi (£) =h ' (/3) *j3 2 , we can conclude there exists a unique 
/3 0 >0 such that h'(Po)=0 and h(j3) attains a minimum at Po; 
and h(/ 3 ) is a monotone increasing function of p when P>Po, 
and h(P) is a monotone decreasing function of p when P<Po- 

Transferring the result from h(/3) back to T(/3) , we conclude 
the Theorem 2 . 


3. Proof of Theorem 3: 

From ( 6 ) , we have 

R c(0) = EXP(-c(r) T ^~|^~ tQ) - ) 

Using the notations defined in the Notation Section, we have 



(3.1) 


Using the first order Taylor expansion, R.H.S. (right hand side) 


of (3.1) = 


[tn/tm). ( (T-etp) /tm)^" 1 


, where O<0<1 

: P (tn/t m ) (T/t-m) ff~l _ (3.2) 


... _ . , -ijii i r.r. \ P l ) ^ H i wn ; \ / ’-m/ = rr 

Considering P>1, we get — c(r) “ SziP z 

Define zi(/3) = Ln (z (/3) tm/to) — Ln/?+( 0 -l) Ln(T/tm) -Ln(Ezi$) . 
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The first derivative and the second derivative of z±((3) are 


, /m 1 . T ,T . Ezi^Lnzi 
zi’(^) = J5 + Ln(^) - 

zi”(/3) = - 1,- (Szi^CLnzQM fSzjff) -(Sz^Lnzj) 


(3.3) 

(3.4) 


By the Cauchy Inequality, the second term of R.H.S. of (3.4) < 0. 
So zi" ()3) <0. Therefore zi'(/J) is a strictly monotone 
decreasing function of /? . Because of the given conditions, we also 
have 

- l + Ln(^) - > 0 

and ^1 j^zi' (/3) = 0 + Ln(T/t m ) - 0 < 0. 

So there exists a unique (3q>0 such that zi'(/3)=0 and zi(/3) 
attains a maximum at /3 o - But noticing 

z(/J) = (t 0 /tm)EXP(zi(/3)), 

we could extend the result to z(j3) . So we have 


-Ln(Rc (/?).)< gtfl Tg-1 (fll < _, fl . = ^QtoT^ 
c ( r ) - Sti^ ~ Z(/3) - Z(/?o) " st .^o 


,/3o"l 


for all /3>1. 


J30-1 


That is, R c (/3) > EXP(-c(r) - Q ^ - Q - ^ ), for all /3>1. 


St- 


P 0 


4. Proof of Theorem 4: 

From (6) , we have 

-LnRr. Tr-ff ~ (Tr-tn)fl = (TrVt m ) ff- f (T^-t n ) /tm) £ f4 ^ 

C(T) Stjpf “ EzjP 

First, we have to show, for any R c , c(r), to, j3 , tj.'s and 

under the given conditions (8), a unique T c =T c (/3) can be solved for 
from (4.1) . 

Define f(T c ) = T C P- (T c -t 0 ) P-XtiP (=j^) . 

Taking the derivative with respect to T c , we have 
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f ' (T c ) =(3T c P-l-p (T c -to) P-^O . 

We also have 1 (T c ) = +® >0 and 

li^j f(T c ) = t 0 ff - Stiff ( -^ ) < 0, because of (8). 

Therefore, a unique solution T c =T c (ff) exists for the equation (4.1). 

Using a Taylor expansion, we have 

-LnRc _ ff(to/t -m) ( (T c -9to) /t m ) ff — 1 where 0<G<1 (4.2) 

c ( r ) Sziff 


Solving for Tc(ff) from (4.2), we have 


T^(ff) _ , -LnRr. Szjff tm. l/(ff-l) , 
t m " ( c(r) ff t 0 ) + 


m 


> , ~ LnRc Sz 'i /3 £nhV(ff- x )= hi ( B) 
£ 1 c ( r ) p to 1 - w 


Let W = t m (-LnR c )/c(r)/to, (constant, independent of ff) and 

LnW+Ln(Zziff) -Ln6 
wi(ff) ■ Lnw ( B ) = ^1 


The first derivative of wi(6) is 

-LnW S |^fl nZi (ff~ 1 )- Ln (S z iff) 

Wl '(ff) = (/3-1) 2 + ^ (ff-1) 2 


( 6-1) / 6~Ln6 
(ff-l) : 


-l+l//3-LnW+Ln/3+ (6~1) -Ln(Zziff ) 

= ( ff - 1) 2 

Let the numerator of R.H.S of (4.3) = W2 (ff) • Then 

6-1 . SziffLnzi ... (Zz-j ff (Lnz-j) 2 ) (Szi ^-(ZzjffLnzjUL 2 

w 2 ’(ff) = yr + sziff + (/ 3 - 1 ) (Sziff ) 2 - 

Zz-iffLnzi 
Zz iff 


(4.3) 


ff-1 . , fl - . f Zz i ff (Lnz i ) 2 ) (Zzjff)-(ZziffLnzi) 2 
B 2 X) (Zziff) 1 


(4.4) 


Using the Cauchy Inequality and noticing ff-l>0, we have W2 , (ff) > 0 
So W 2 (ff ) is a strictly monotone increasing function of 6- 

Using the given condition: Rq < EXP( -c ( r) . ), we have 
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and 


W 2 (jS) = -LnW - Ln(Ezi) < 0 
pl^oo w 2(0) = + 00 > 0. 

So there exists a unique (3q>1 such that W2(jSo) = 0 an< ^ 

W2(£) < 0 when (3 < (3q; W 2 (£) > 0 when (3 > Po- 

But wi 1 ((3) , so the same conclusions hold for wi'(/?) for /3>1. 

Therefore, wi(/?) attains a minimum at /3o* Hence w(£)=EXP(wi ((3) ) 
attains a minimum at /3 0 - That is, 

^ > w(/ 3 ) > .(ffO)-(=MB ic £2i^ &J, Vtfo-U . 

Therefore we have: T c (/S) ^ ( -Ln (Rc) — ( r^o^O ^ / for a H /3>1. 
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